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Section  I 

Introduction  and  Summary 

In  this  report,  we  examine  the  effect  of  spall  on  regional  and  teleseismic 
waveforms  using  two-dimensional  nonlinear  numerical  simulations  of  explo¬ 
sions.  Each  of  the  following  sections  covers  a  different  topic  related  to  this 
theme.  Sections  II  and  III  are  summaries  of  material  that  will  be  described  in 
more  detail  in  scientific  reports  which  are  now  in  preparation. 

In  Section  II,  we  develop  a  simple,  physical  model  for  spall  by  comparing 
the  complete  two-dimensional  waveform  with  the  waveform  generated  by  a 
one-dimensional  explosion  plus  a  shallow  tension  crack  in  the  same  layered 
medium.  This  is  an  extension  to  higher  frequencies  of  the  spall  model  derived 
by  Day,  et  al.  (1983).  The  waveform  for  a  one-dimensional  explosion  is  sub¬ 
tracted  from  the  waveform  for  the  two-dimensional  simulation,  and  the  residual, 
which  contains  all  of  the  complex  free-surface  interactions,  is  interpreted  as  the 
waveform  generated  by  the  tension  crack.  The  tension  crack  is  parameterized 
by  its  radius,  depth,  and  a  distribution  of  takeoff  velocities  over  its  surface. 

We  find  that  the  P-waves  generated  by  the  two-dimensional  simulation  are 
modeled  very  well  by  the  simple  tension  crack  plus  explosion  model.  The 
resulting  parameters  for  the  spall  model  are  consistent  with  the  limited  set  of 
near-field  obsen/ations  of  spall.  The  peak  spall  velocity  is  found  to  be  much 
higher  for  underburied  explosions  than  for  normal  or  overburied  explosions  The 
spall  depth  decreases,  but  only  slowly,  as  the  explosion  depth  decreases.  The 
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model  underestimates  the  shear  waves  generated  in  the  two-dimensional  calcu¬ 
lations,  indicating  that  the  source  of  shear  waves  is  more  complex  than  the  sim¬ 
ple  explosion  plus  tension  crack  model. 

The  spall  model  is  found  to  have  a  narrowband  spectrum  proportional  to  co^ 
at  low  frequencies  and  at  high  frequencies.  The  spall  spectrum  is  peaked 
at  a  frequency  inversely  proportional  to  the  spall  dwell  time.  A  typical  range  for 
the  peak  frequency  is  from  0.5  to  5  Hz,  which  puts  it  in  the  middle  of  the  fre¬ 
quency  range  where  short  period  magnitudes  are  commonly  measured. 

In  Section  III,  we  use  the  model  developed  in  Section  II  to  generate 
regional  seismograms  for  an  explosion  plus  spall  source.  The  source  function 
derived  from  the  spall  model  varies  as  a  function  of  takeoff  angle,  or  equivalent 
phase  velocity,  so  the  effective  source  functions  are  different  for  each  regional 
phase.  In  particular,  the  lower  the  phase  velocity,  the  more  narrowband  the 
spall  source,  so  the  Lg  generated  by  the  spall  source  is  significantly  more  nar¬ 
rowband  than  the  Pg,  Pn,  or  teleseismic  P-waves. 

In  high  velocity  structures,  such  as  the  model  we  used  for  the  simulation  of 
Shagan  River  area  explosions,  a  spherically  symmetric  explosion  source  is  a 
very  poor  generator  of  Lg.  In  contrast,  the  tension  crack  is  a  good  generator  of 
Lg,  so  the  Lg  generated  by  spall  completely  dominates  the  regional  Lg  syn¬ 
thetic.  For  a  125  kiloton  simulation,  we  estimate  an  mt,{Lg)  of  6.2  for  the  spall 
source  and  only  4.9  for  the  explosion  source.  The  regional  phases  Pn  and  Pg, 
on  the  other  hand,  are  comparable  in  amplitude  for  the  spall  and  explosion 
sources.  This  is  expected  since  the  spall  signal  has  the  effect  of  partially 
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canceling  the  elastic  pP  signal  and  replacing  it  with  a  delayed  and  lower  fre¬ 
quency  phase. 

In  Section  IV,  we  examine  the  effect  of  variations  in  earth  structure  on  the 
generation  of  Lg  by  explosion  and  spall  sources.  This  study  was  motivated  by 
the  results  of  Section  III;  in  particular  by  the  need  to  assess  the  robustness  of 
the  rather  surprising  conclusion  that  nearly  all  of  the  Lg  synthetic  for  an  explo¬ 
sion  is  due  to  spall.  This  result  was  tested  by  generating  synthetic  regional 
seismograms  in  a  set  of  earth  models  derived  by  several  authors  for  Eastern 
Kazakhstan  near  the  Shagan  River  test  site.  The  result  was  found  to  be 
independent  of  any  of  these  structures.  The  characteristic  of  these  structures 
that  is  responsible  for  this  effect  is  a  high  source  region  P-wave  speed,  in 
excess  of  the  Lg  phase  velocity.  Under  this  condition,  explosion  Lg  can  be 
thought  of  as  being  generated  by  S*,  a  nongeometric  arrival  caused  by  the  cur¬ 
vature  of  the  wavefront.  Since  this  phase  decays  exponentially  with  source 
depth,  it  is  a  very  weak  generator  of  Lg,  and  the  resulting  Lg  will  be  swamped 
by  the  Lg  from  spall  or  any  other  deviatoric  part  of  the  source. 

In  contrast,  for  a  structure  with  low  velocities  at  the  source,  Lg  generation 
by  an  explosion  source  is  comparable  to  the  generation  of  Lg  by  the  spall 
source.  The  Lg  phase  generated  by  the  spall  or  CLVD  source  is  more  complex 
than  the  Lg  phase  generated  by  the  point  explosion  source:  however,  the  peak 
amplitudes  are  similar,  and  Lg/Pg  peak  amplitude  ratios  are  approximately 
equal. 
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It  can  be  concluded  from  the  results  of  Sections  III  and  IV  that  a  point 
explosion  source  is  an  unrealistic  model  for  the  generation  of  Lg  for  high  velo¬ 
city  structures.  Some  physical  effect  must  break  the  spherical  symmetry  of  the 
source  in  order  to  generate  the  observed  Lg  amplitudes.  Spall  is  one  mechan¬ 
ism  that  breaks  the  spherical  symmetry  at  a  level  sufficient  to  match  the  data. 
Other  mechanisms,  such  as  near  source  scattering,  may  also  generate  enough 
Lg  to  explain  the  data  (see  Section  VI). 

In  Section  V,  we  examine  the  regional  seismic  signals  of  the  NTS  explo¬ 
sions  Duryea  and  Buteo.  The  significance  of  these  explosions  is  that  they  are 
at  the  same  location,  but  with  very  different  scaled  depths  of  burial.  Duryea 
was  at  normal  scaled  depth,  while  Buteo  was  significantly  overburied.  Since 
spall  should  be  greater  for  a  shallow  source  than  for  a  deep  source,  a  com- 
pahson  of  these  two  explosions  provides  an  opportunity  to  see  if  the  generation 
of  Lg  is  related  to  the  amount  of  spall.  Analysis  of  the  data  shows  that  there  is 
no  measurable  differerice  between  the  regional  signals  from  these  two  signals 
other  than  that  due  to  standard  yield  scaling,  so  the  difference  in  spall  does  not 
appear  to  have  resulted  in  a  difference  in  Lg.  Unfortunately,  since  these  explo¬ 
sions  took  place  in  a  low  velocity  structure,  where  Lg/Pg  ratios  are  expected  to 
be  similar  for  explosion  and  spall  sources,  it  is  not  possible  to  separate  the  spall 
source  from  the  explosion  source  in  the  data.  A  similar  comparison  in  a  high 
velocity  structure  would  be  definitive;  however,  no  opportunity  for  such  a  com¬ 
parison  exists  at  present. 
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In  Section  VI,  we  summarize  some  of  the  additional  research  that  is  in  pro¬ 
gress  under  this  contract.  One  important  area  of  research  is  the  implementa¬ 
tion  of  damage  mechanics,  as  developed  by  Sammis  and  Ashby,  in  S-CUBED 
nonlinear  finite  difference  codes.  The  significance  of  this  damage  model  is  that 
it  is  scale  dependent,  and  may  be  able  to  explain  observed  strength  differences 
between  laboratory  and  in  situ  rock.  A  second  area  of  research  is  the  effect  of 
scattering  on  Lg.  To  this  end,  we  are  performing  a  finite  difference  simulation 
of  the  Salmon  experiment,  using  a  realistic  model  of  the  highly  nonuniform 
structure  near  this  explosion. 
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Section  II 

A  Simple  Physical  Model  for  Spall  from  Nuclear  Explosions  Based  Upon 
Two-Dimensional  Nonlinear  Numerical  Simulations 

Introduction 

Several  models  have  been  proposed  tor  the  spall  process  that  accom¬ 
panies  most  nuclear  explosions.  The  models  are  very  difficult  to  validate 
directly  because  (1)  the  features  of  the  process  that  can  be  measured,  such  as 
vertical  motion  at  the  surface,  are  observed  at  only  a  few  locations  above  the 
shot,  and  (2)  there  are  few  observations  of  features  such  as  the  depth  and  area 
of  the  spall  zone  and  the  distribution  of  motions  within  the  zone.  Two- 
dimensional  nonlinear  simulations  of  the  explosion,  which  include  the  physics  of 
the  free-surface  interactions,  provide  an  opportunity  to  examine  spall,  and  to 
find  simple  models  to  represent  it.  In  this  section,  we  describe  the  work  by 
Barker  and  Day  (1990),  who  present  a  li  ’ear  model  for  the  spall  process  and 
find  its  parameters  by  comparing  with  the  seismic  waves  from  explosion  simula¬ 
tions. 

Using  a  form  of  the  elastodynamic  representation  theorem,  we  have  calcu¬ 
lated  the  short-period  body  waves  emanating  from  the  source  zone  of  the 
numerical  simulations.  The  implications  for  teleseismic  magnitude  measure¬ 
ments  were  discussed  in  Day,  et  al.  (1986)  and  McLaughlin,  et  al.  (1988).  By 
comparing  these  P  and  SV  wavefields  with  one-dimensional  nonlinear  calcula¬ 
tions  (which  do  not  include  the  nonlinear  effects  of  the  free-surface),  we  can 
isolate  that  part  of  the  wave  field  due  to  the  nonlinear  interaction  with  the  free- 
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surface.  We  find  that  a  simple  tension  crack  with  an  opening  that  propagates 
with  the  pP  arrival  from  the  explosion  fits  P-waves  from  simulations  of  both 
Pahute  Mesa  and  Shagan  River  tests.  However,  the  model  which  fits  the  P- 
wave  radiation  generates  SV-waves  that  are  too  small,  which  indicates  that 
there  are  additional  sources  of  shear  waves  in  the  finite  difference  simulations. 

The  tension  crack  model,  proposed  in  its  original  form  by  Day,  et  al. 
(1983),  is  a  physical  model  whose  parameters  can  be  compared  directly  to  field 
observations.  The  parameters  of  the  model  are  spall  depth  and  area,  detach¬ 
ment  velocity,  and  momentum.  The  parameters  inferred  from  the  simulations 
are  consistent  with  published  estimates  of  these  values.  The  time  dependence 
of  the  far-field  waveforms  is  a  natural  consequence  of  the  model  and  requires 
no  ad  hoc  choice  of  time  history  to  include  the  effects  of  source  finiteness,  the 
importance  of  which  has  been  pointed  out  by  Stump  (1985).  Although  there  is 
a  trade-off  between  the  parameters  of  the  model,  the  amplitudes  and 
waveforms  from  the  two-dimensional  (2-D)  simulations  tightly  constrain  the  set 
of  parameters  which  fit  the  simulations. 

The  form  for  the  linear  source  representation  of  the  spall  model  is  such 
that  it  can  easily  be  added  to  a  one-dimensional  (1-D)  explosion  source  to  com¬ 
pute  regional  and  teleseismic  synthetic  seismograms.  In  a  companion  paper, 
McLaughlin,  et  al.  (1990),  use  the  model  to  compute  the  effects  of  spall  on 
synthetic  regional  explosion  seismograms.  It  is  found  that  the  spall  contribution 
to  the  Lg  signal  is  comparable  to  the  direct  explosion  contribution.  Their  work  is 
summarized  in  Section  III  of  this  report. 
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The  Tension  Crack  Model 

We  have  extended  the  circular  tension  crack  model  of  Day,  et  at.  (1983)  to 
include  a  time  dependence  which  describes  the  detachment  and  slap-down  of 
the  spall  volume.  The  spall  volume  in  the  model  is  defined  as  the  cylinder 
whose  radius  is  that  of  the  tension  crack  and  which  extends  from  the  crack  to 
the  surface  (Figure  2.1).  Since  we  wish  to  compare  the  spall  model  with  2-D 
simulations  which  were  done  in  a  layered  medium,  we  developed  the  formalism 
for  computing  the  far-field  body  waves  emanating  from  a  layered  medium.  We 
briefly  describe  the  features  of  the  model  and  refer  the  reader  to  Barker  and 
Day  (1990)  for  derivations  and  details.  The  parameters  of  the  model  are  the 
crack  radius  and  depth,  and  the  distribution  of  detachment  velocities  over  the 
crack.  In  the  complete  model,  as  described  in  Barker  and  Day,  the  tension 
crack  opening  is  tied  to  the  arrival  time  along  the  spall  surface  of  the  reflected 
tensile  wave  pP.  In  addition,  the  initial  opening  velocity  varies  smoothly  over 
the  spall  surface.  The  spall  volume  then  decelerates  under  gravity  until  it 
closes. 

The  general  behavior  of  the  model  can  be  illustrated  for  the  analytically 
simpler  case  in  which  the  crack  opens  simultaneously  over  the  whole  spall  sur¬ 
face,  and  the  initial  velocity  is  uniformly  distributed  over  the  surface  and  over 
the  velocity  range  to  V2.  In  that  case,  we  can  approximate  the  far-field  P- 
wave  displacement  along  the  ray  by 
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^ray(0  ~ 


p^o?R 


[(ap)2-f2]i/2  *  Soit)  ,  -ap  <t  <ap 


(1) 


where  p  is  the  ray  parameter,  a  is  the  radius  of  the  tension  crack,  is  the 
depth  of  the  crack,  a  is  the  P-wave  speed  at  the  source,  is  the  geometric 

spreading  factor,  and  *  denotes  convolution.  The  function  So(0  is  the  time  His¬ 
tory  of  the  mean  acceleration  of  the  spall  volume,  and  is  given  by 


So{t)  =  vm-9^H(0.t,)  + 


■  g- 

V2-V1 


i-V2+ j9t)  AH {t„t2)  ■ 


(2) 


Here,  and  V2  are  the  the  minimum  and  maximum  detachment  velocities  on 

2v, 

the  crack,  the  times,  //= — ,  are  the  corresponding  flight  times,  V  is  the  mean 

Q 

detachment  velocity,  and 


is  a  "boxcar"  turning  on  at  and  off  at  f2-  Equation  (2)  is  shown  schematically 
in  Figure  2.2,  The  displacement  time  function  u^y{t)  is  formed  by  convolving 
Soit)  with  (a^p^  -  ,  which  results  from  the  radial  finiteness  of  the  spall. 

We  note  that  typical  teleseismic  values  of  p  are  less  than  0.1  sec/km,  and 
values  of  a  are  around  1  km,  so  ap  <  0.1  sec,  which  allows  (1)  to  be  further 
approximated  by 


(-va+lsf)  AHM)  I 

V2-V^  4 
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where  nrispaii  is  the  mass  of  the  spall  volume,  and  X  and  |i  are  the  Lame’ 
parameters  at  the  source.  The  inital  5  function  term  in  (4)  is  proportional  to 
rr^spauV  .  the  mean  spall  momentum.  The  dwell  term  is  proportional  to  rVspang, 
the  gravitational  restoring  force.  The  duration  of  the  dwell  is  and  that  of  the 
slap  down  is 

The  spectrum  of  the  tension  crack  is  shown  schematically  in  Figure  2.3.  It 
is  peaked  near  the  frequency  =  Typical  maximum  detachment 

velocities  are  between  1  and  10  m/sec,  so  values  of  t2  are  .2  to  2  sec.  Assum¬ 
ing  the  minimum  velocity  is  nearly  zero,  fpi^  would  be  in  the  range  0.5  to  5. 
Hz,  which  puts  the  peak  right  in  the  teleseismic  and  regional  frequency  bands. 
The  rate  of  roll-off  at  high  frequencies  depends  on  the  spatial  dependence  of 
detachment  velocity,  but  for  any  distribution  which  goes  to  zero  at  the  edge  of 
the  crack,  the  rate  is  at  least  as  fast  as  The  spectrum  in  Figure  2.3  rolls 
off  at  both  high  and  low  frequencies  faster  than  explosion  source  functions,  so 
that  the  effects  of  spall  on  the  seismic  signal  are  limited  to  a  narrow  frequency 
band. 

Comparison  with  Numericai  Simulations 

The  objective  of  this  study  was  to  derive  a  simple,  linear  physical  model 
which  represents  the  spall  process.  Our  strategy  has  been  to  isolate  surface 
interaction  effects  by  comparing  the  wavefields  from  2-D  calculations,  U2_o. 
with  those  from  1-D  calculations  done,  u^_q,  in  the  same  source  materials  (but 
in  an  inifinite  medium).  Our  working  assumption  is  that  the  difference  in  the 
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Log  Frequency 


Figure  2.3.  Spectrum  of  far-field  P-waves  from  the  tension  crack  model. 
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wavefields  represents  the  free-surface  interaction  effects.  In  fact,  the  procedure 
has  been  to  find  the  difference  in  the  displacements  U2-D~^^-D< 
parameters  of  the  spall  model  (described  in  the  previous  section)  which  have 
displacements  Uspan  which  best  satisfy  Uspan  =  U2-D~<^^-D■  Alternatively,  we 
could  have  found  Uspaii  which  satisfies  Ui_o+Uspa//  =  ^2-0-  We  found  the  first 
approach  to  be  more  direct. 

To  find  the  ground  motions  from  the  2-D  calculations,  we  used  the 
methods  described  in  Rodi,  et  al.  (1978),  Bache,  et  al.  (1982),  and  Day, 
et  al.,  (1983,  1986).  The  2-D  simulations  themselves  are  described  in  Day 
et  al.  (1986).  The  elastic  properties  for  the  Pahute  Mesa  and  Shagan  River 
simulations  are  shown  in  Table  2.1.  The  method  for  computing  the  body  waves 
entails  evaluating  a  form  of  the  elastodynamic  representation  theorem,  which 
gives  the  motions  in  terms  of  spatial  and  temporal  convolutions  of  displace¬ 
ments  and  stresses  monitored  on  a  surface  surrounding  the  nonlinear  zone  with 
Green’s  functions  and  their  spatial  gradients. 
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Table  2.1 

Elastic  Structure  for  the  Shaqan  River  Simulations 

Layer 

Compressional 

Shear 

Density 

Thickness 

Velocity 

Velocity 

(m) 

(m/sec) 

(m/sec) 

Kg/ 

oo 

5018 

2789 

2700 

Elastic  Structure  for  the  Pahute  Mesa  Simulations 

Compressional 

Shear 

Density 

Velocity 

Velocity 

■Bn 

(m/sec) 

(m/sec) 

Kg/ 

112.5 

1208 

661.4 

1600 

457.5 

2025 

1109 

1950 

OO 

2887 

1581 

2000 

Shagan  River  Simulations 

We  begin  with  the  Shagan  River  simulations,  which  were  all  done  at  a  yield 
of  125  KT  and  at  four  depths  of  burial  (DOB):  an  over  buried  depth  (980  m),  an 
optimally  buried  depth  (680  m),  a  depth  just  below  cratering  (no  ejecta)  (300  m), 
and  a  depth  which  causes  cratering  (200  m).  The  scaled  depths  of  burial  are 
196,  120,  60  and  40  m  The  results  for  P-waves  were  presented  in 

Day,  et  al.  (1986),  and  we  include  them  here  for  purposes  of  comparing  with 
the  spall  model.  As  can  be  Seen  in  Table  2.1,  the  earth  model  is  a  half-space, 
which  reflects  our  lack  of  detailed  knowledge  of  the  structure  rather  than  making 
a  statement  of  homogeneity  at  the  site.  The  far-field  P  and  SV  displacements 
are  shown  for  a  take-off  angle  of  10°  in  Figures  2.4  and  2.5.  In  our  detailed 
report  (Barker  and  Day,  1 990),  we  make  comparisons  at  three  take-off  angles, 
10°,  20°  and  30°,  which  are  representative  of  teleseismic,  intermediate  and 


15 


Id,  dob-300 


Comparison  of  the  far-field  SV-waves  from  the  1-D  (solid  lines)  and  2-D  (dashed  lines)  Shagan 
River  simulations  for  four  depths  of  burial  at  a  take-off  angle  of  10” . 
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regional  slownesses  (phase  velocities)  for  this  source  structure.  These  com¬ 
parisons  ensure  that  the  spatial  dependence  is  modeled  correctly.  In  Figures 
2.4  and  2.5,  the  motions  from  the  2-D  simulations  are  overlain  with  1-D  simula¬ 
tions  in  which  the  nonlinear  properties  at  the  source  are  the  same  as  the  2-D 
values.  The  1-D  source  (RDP)  was  extracted,  and  the  body  waves  were  calcu¬ 
lated  for  a  linear  elastic  medium.  Thus,  for  the  Shagan  River  half-space  model, 
the  1-D  waveforms  have  the  direct  P  phase,  as  well  as  the  pP  and  pS  elastic 
reflections.  These  signals  correspond  to  down-going,  teleseismic  ray  paths,  but 
do  not  include  the  effects  of  mantle  or  crustal  propagation,  anelastic  attenua¬ 
tion,  or  recording  instrumentation. 

The  rise  time  and  amplitude  of  the  first  peak  in  the  P-waves  (Figure  2.4)  is 
very  nearly  the  same  for  the  1-D  and  2-D  calculations.  For  the  shallow  calcula¬ 
tions  (200  and  300  m)  the  1  -D  and  2*D  waveforms  diverge  after  the  peak  and 
are  quite  different  at  later  times,  indicating  that  the  cratering  and  near-cratering 
processes  are  indeed  different  from  elastic  pP  reflections.  For  the  deeper 
DOB’S,  680  and  980  meters,  the  apparent  pP  phase  on  the  2-D  records 
appears  to  be  smaller,  to  have  longer  period,  and  to  be  delayed  relative  to  the 
elastic  pP  (1-D)  case. 

The  SV-waves  (Figure  2.5)  show  much  greater  differences  between  the  1- 
D  and  2-D  cases.  The  2-D  solutions  have  a  direct  S-wave,  due  to  vertical 
asymmetries  in  the  source,  which  is  not  in  the  1-D  solutions.  The  main  peaks 
in  the  2-D  SV-waves  are  larger  than  the  1-D  peaks  at  all  four  DOS’s,  and  are 
about  twice  as  big  for  the  300,  680  and  980  depths.  The  duration  of  the  main 
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peaks  is  also  greater  for  the  2-D  cases.  A  large  negative  swing  occurs  on  the 
2-D  waveforms,  which  is  much  smaller  on  the  1-D  signals. 

As  discussed  above,  we  computed  the  2D-1D  difference  waveforms  and 
matched  them  to  those  from  the  spall  model.  The  parameters  of  the  spall 
model  which  most  closely  matched  the  numerical  simulations  are; 


Table  2.2 

Spall  Model  Parameters  for  the  Shagan  River  Simulations 

Depth 

Crack 

Crack 

Minimum 

Maximum 

of  Burial 

Depth 

Radius 

Detachment 

Velocity 

Detachment 

Velocity 

•^exp 

a 

V'l 

V2 

(m) 

_ Im) _ 

(m) 

(m/sec) 

(m/sec) 

600 

1.1 

20.0 

150 

600 

1.1 

15.0 

680 

600 

1.1 

4.0 

980 

200 

600 

1.1 

1.5 

In  Figure  2.6,  we  compare  the  P-wave  spall  model  waveforms  with  the  2D-1 D 
waveforms  for  each  of  the  DOB’s.  In  general,  the  comparisons  are  quite  good, 
both  in  amplitude  and  shape.  The  parameter  that  varies  the  most  in  Table  2.2 
is  the  maximum  detachment  velocity  V2  .  Recall  that  in  Equation  (4),  V2  enters 
the  solution  in  two  ways.  First,  the  5  function  term,  which  causes  the  first  peak 
in  the  waveforms,  is  proportional  to  the  mean  detachment  velocity  7.  Second, 
the  duration  of  the  sig''al  (end  of  slap-down)  is  time  t2  =  2v2/g.  As  the  DOB 
decreases,  the  2D-1D  difference  waveforms  increase  in  amplitude  and  duration. 
Thus,  as  the  DOB  decreases,  V2  is  required  to  increase.  We  note  that,  even  in 
the  case  where  cratering  occured  (200  m),  the  tension  crack  model  provided  a 
good  representation  of  the  P  radiation.  The  shapes  of  the  signals  for  the  980 


t9 


SSS-TR-90-11536 


meters  case  agree  well.  For  the  680  meters  case,  the  agreement  is  good 
except  for  an  additional  inflection  in  the  later  parts  of  the  2-D  waveforms  that  is 
not  modeled  by  the  tension  crack. 

The  peak  detachment  velocity  is  inferred  from  field  observations  from 
ground  motion  records  at  or  near  ground  zero.  We  can  do  an  analogous  meas¬ 
urement  by  examining  the  vertical  velocity  in  the  finite  difference  simulation  at 
the  free-surface  of  the  grid.  We  find  that  the  values  of  V2  in  Table  2.2  agree 
closely  with  the  ground  zero  velocities  in  the  simulations. 

Using  the  parameters  in  Table  2.2,  we  computed  the  SV-waves  from  the 
tension  crack  model  and  compared  them  with  the  corresponding  2D-1D 
difference  SV-waves.  As  can  be  seen  in  Figure  2.7,  the  comparisons  are  not 
as  good  as  the  P-waves,  especially  at  the  deeper  DOS’s,  where  the  amplitudes 
predicted  by  the  model  are  too  small  by  factors  of  two  to  three.  Attempts  to  find 
a  set  of  parameters  which  would  Improve  the  SV  comparisons  without  degrad¬ 
ing  the  P  comparisons  were  unsuccesful.  We  note  that  in  all  cases,  the  tension 
crack  models  derived  from  the  P-waves  account  for  only  a  fraction  of  the  SV 
radiation  from  the  nonlinear  simulations. 

The  results  of  the  SV  comparisons  suggest  that  additional  sources  of  SV 
radiation  be  incorporated.  We  therefore  tried  some  additions  to  the  model  to 
improve  this  situation.  First,  we  hypothesized  that  processes  transpired  in  the 
simulations  that  could  be  represented  as  modifications  to  the  isotropic  part  of 
the  moment  tensor.  We  chose  to  represent  this  as  a  compensated  linear  vector 
dipole,  CLVD  (  e.g.,  Knopoff  and  Randall,  1970).  The  CLVD  represents  a 
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process  akin  to  squeezing  a  vertical  tube  of  toothpaste.  We  also  hypothesized 
that  additional  sources  of  shear  waves  may  be  emanating  from  vertical  cracks 
which  are  the  boundaries  of  a  cylinder  which  moves  above  the  horizontal  ten¬ 
sion  crack.  That  is,  we  imagine  that  the  material  above  the  tension  crack 
moves  up  and  down  as  a  unit,  shaped  as  a  cylinder  whose  base  is  the  tension 
crack.  The  relative  motion  along  the  vertical  sides  would  act  as  a  distributed 
shear  dislocation.  In  Barker  and  Day  (1990),  we  compare  the  SV-waves  from 
these  two  models  with  the  2D-1D  waveforms.  The  CLVD  source  when  added 
to  the  tension  crack  source  in  various  proportions  produced  SV-waves  that  were 
large  enough  to  fit  the  2D-1D  difference  waveforms,  and  the  shape  of  the 
waveforms  agreed  well.  Since  the  CLVD  source  generated  small  amounts  of 
P-waves,  the  goodness  of  fit  to  the  P-waves  seen  above  was  not  affected. 
However,  the  proportions  were  different  for  each  DOB,  which  effectively  added 
another  free  parameter  whose  interpretation  was  ambiguous.  The  distributed 
shear  dislocation  source  (the  moving  cylinder  described  above)  produced  SV- 
waves  that  were  too  small,  and  so  it  appears  that  the  cylindrical  source  is  not  a 
likely  candidate  for  additional  SV  radiation. 

Pahute  Mesa  Simulations 

The  2-D  and  1-D  P-wave  signals  for  the  Pahute  Mesa  simulations  are 
overlain  in  Figure  2.8  for  a  take-off  angle  3.3°.  The  calculations  were  done  at 
DOB’S  of  200,  680  and  980  m  for  the  Pahute  Mesa  model,  but  not  at  300  m. 
The  first  peaks  of  the  1-D  and  2-D  signals  align  closely.  For  DOS’s  680  and 
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2.8.  Comparison  of  the  far-field  P-waves  solutions  from  the  1-D  (solid  lines)  and  2-D  (dashed  lines) 
Pahute  Mesa  simulations  for  three  depths  of  burial  at  a  take-off  angle  of  3.3* . 
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980,  the  apparent  pP  from  the  2-D  waveforms  is  lagged  about  0.2  seconds  rela¬ 
tive  to  the  elastic  pP  in  the  1-D  waveforms. 

As  with  the  Shagan  River  simulations,  we  fit  the  tension  crack  model  to  the 
2D-1D  difference  time  series.  The  results  are  shown  in  the  comparison  plot  in 
Figure  2.9.  The  fits  of  the  model  to  the  differential  2D-1D  signals  are  generally 
good  for  the  three  DOB'S  and  take-off  angles.  The  best  fitting  parameters  of 
the  model  are 


Table  2.3 

Spall  Model  Parameters  for  the  Pahute  Mesa  Simulations 

Depth 

Crack 

Crack 

Minimum 

Maximum 

of  Burial 

Depth 

Radius 

Detachment 

Velocity 

Detachment 

Velocity 

^exp 

a 

‘'i 

V2 

(m) 

(m) 

(m) 

(m/sec) 

(m/sec) 

200 

150 

400 

1.1 

22 

680 

200 

2200 

1.1 

2 

980 

150 

2200 

1.1 

2 

There  are  several  differences  between  the  model  parameters  derived  here 
and  those  for  the  Shagan  River  simulations.  First,  we  found  for  the  Pahute 
Mesa  runs  that  the  detachment  velocity  must  be  distributed  over  the  crack  with 
the  largest  values  near  the  center  in  order  to  match  the  change  in  2D-1D  sig¬ 
nals  with  take-off  angle.  This  is  in  contrast  to  the  Shagan  River  case,  where  the 
solutions  were  insensitive  to  the  form  of  the  distribution.  This  difference  is  due 
to  the  velocity  structures,  where  the  Shagan  River  wave  speeds,  and  hence  the 
phase  velocities,  are  about  twice  those  of  Pahute  Mesa. 
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Another  feature  of  the  Pahute  Mesa  simulations  that  differs  from  the 
Shagan  River  runs  is  the  large  change  in  crack  radius  between  the  case  for 
DOB=200  and  the  other  two  DOS’s.  This  is  presumably  due  to  the  shallow 
layering  in  the  Pahute  Mesa  earth  model.  Although  the  tension  crack  had  a 
maximum  detachment  velocity  (22  m/sec)  similar  to  Shagan,  it  occurred  over  a 
much  smaller  radius  (400  m). 

The  results  of  the  comparisons  of  the  SV-waves  from  the  tension  crack 
model  in  Table  2.3  are  similar  to  those  for  the  Shagan  River  simulations  dis¬ 
cussed  above.  The  tension  crack  generates  SV-waves  that  are  too  small. 

Comparisons  with  Observations 

From  observations  of  surface  ground  motions  and  of  physical  manifesta¬ 
tions  of  Pahute  Mesa  tests,  several  authors  have  estimated  the  parameters  of 
the  spall  process.  In  Table  2.4,  we  compare  their  results  with  those  for  the  ten¬ 
sion  crack  model  at  the  optimal  DOB  (680m). 
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Table  2.4.  Spall  parameters  from  field  observations  and 
from  the  tension  crack  model  for  Pahute  Mesa 


Parameter 

Tension 

Crack 

m 

Viecelli 

(1973) 

Stump 

(1985) 

Rawson 

(1988) 

Maximum 

Velocity 

(m/sec) 

2 

m 

8 

2200 

1750 

500 

500- 

2665 

IWTHrlSB 

Hi 

9.2 

3.5 

0.58 

1 

1 

5.3 

3.5 

1.2 

H 

H 

200 

110 

100- 

400 

The  authors  cited  in  the  table  typically  expressed  their  results  as  scaled  values, 
in  which  case  we  used  125  KT,  the  yield  in  the  simulations.  In  addition,  the 
values  in  Table  2.4  are  the  average  values  given  by  the  authors.  The  mean 
detachment  velocity  (1.55  m/sec)  was  used  to  compute  the  tension  crack 
momentum,  rather  than  the  mean  velocity.  The  momentum,  spall  radius  and 
depth  from  the  tension  crack  model  lie  within  the  spread  of  those  inferred  from 
observations.  The  mass  is  at  the  high  end  while  the  maximum  velocity  is  at  the 
low  end. 

It  is  very  difficult  to  estimate  the  parameters  in  Table  2.4  because  of  a  pau¬ 
city  of  direct  observations.  There  are  typically  few  surface  ground  motion  sen¬ 
sors,  and  no  sub-surface  observations  of  ground  motion  or  other  phenomena 
such  as  cracking.  Inference  of  the  spall  process  beyond  estimating  apparent 
pP  amplitude  and  travel  time  from  far-field  recordings  has  not  been  successful. 
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With  these  limitations  in  mind,  it  can  be  said  that  the  parameters  of  this  study 
are  consistent  with  published  field  observations. 

Summary  and  Conclusions 

We  show  that  a  simple  model  can  be  constructed  which  generates  the  far- 
field  P-waves  of  a  two-dimensional  nonlinear  calculation  which  includes  the 
effects  of  the  free-surface  for  a  range  of  depths-of-burial.  The  model  which  fits 
the  P-waves  from  the  simulations  underestimates  the  SV-waves.  The  model, 
proposed  in  its  original  form  by  Day,  et  al.  (1983),  is  based  on  using  a  tension 
crack  which  opens  due  to  the  tension  wave  from  the  free-surface.  The  material 
over  the  crack  travels  up  with  the  impulse  of  the  tension  wave  and  returns 
under  the  influence  of  gravity.  Stump  (1985)  modified  the  original  model  to 
include  the  effects  of  source  by  adding  an  empirical  time  function  based  on 
chemical  explosions.  The  modified  model  was  used  by  Taylor  and  Randall 
(1988)  to  model  regional  seismograms.  In  our  formalism  presented  in  this 
report,  we  include  the  effects  of  source  finiteness  and  those  of  crustal  rever¬ 
berations.  A  time  dependence  is  a  natural  consequence. 

The  parameters  of  the  model  compare  favorably  with  observations  based 
on  field  data.  In  addition,  McLaughlin,  et  al.  (1990),  use  the  model  to  compute 
regional  seismograms  and  find  that  in  the  Lg  bandwidth,  the  spall  contribution 
should  be  comparable  to  or  greater  than  that  of  the  explosion  signal  alone. 
Since  the  spall  model  which  fits  the  P-waves  in  our  study  generates  SV-waves 
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that  are  too  small,  the  spall  model  may  lead  to  an  underestimate  of  the  Lg  from 
spall.  Spall,  therefore,  appears  to  be  a  very  significant  source  of  Lg. 
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Section  III 

Implications  of  Explosion  Generated  Spall  Models: 

Regional  Seismic  Signals 

Spall  is  a  frequently  observed  phenomenon  associated  with  the  free  sur¬ 
face  reflection  of  an  explosion  generated  compressional  wave.  The  compres- 
sional  wave  originating  from  an  underground  explosion  becomes  a  tensional 
wave  upon  reflection  from  the  free-surface.  If  the  rock  near  the  free-surface  fails 
under  this  tension  then  the  volume  of  material  above  the  tensile  failure  goes 
into  ballistic  free-fall.  This  free-fall  is  often  recognized  by  acceleration  records 
above  the  explosion  that  exhibit  -1  g  dwells  (free  fall)  followed  by  a  sudden  slap 
down.  During  this  free  fall  the  motion  of  the  material  is  largely  governed  by  the 
force  of  gravity.  A  number  of  authors  have  gathered  data  concerning  spall  from 
both  underground  nuclear  explosions  and  contained  chemical  blasts  (Eisler  and 
Chilton,  1964;  Viecelli,  1973;  Sobel  1978;  Stump,  1985;  Patton,  1990).  The  Day, 
et  at.  (1983)  seismic  representation  for  spall  has  become  a  widely  accepted 
model  for  the  source  representation  of  this  nonlinear  phenomenon.  This  model 
has  been  modified  and  used  by  Stump  (1985),  Patton  (1988),  and  Taylor  and 
Randall  (1989)  to  model  near-field  and  regional  seismic  phases  due  to  spall. 
The  models  for  spall  considered  here  are  based  on  the  original  Day,  et  a!. 
(1983)  tension  crack  representation  for  spall  and  modified  to  include  a  distribu¬ 
tion  of  spall  over  the  zone  of  tensile  failure  as  parameterized  by  Barker  and  Day 
(1990). 
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Assessing  the  importance  of  spall  in  exciting  regional  seismic  signals 
requires  a  physically  reasonable  and  consistent  model  for  the  forces  and/or 
moment  tensor  expansion  of  those  forces.  These  representations  may  then  be 
used  to  make  predictions  of  the  spall  signal  and  to  invert  for  possible  spall 
models  using  far-field  seismic  data. 

Section  II  of  this  report  summarizes  the  results  from  modeling  2-D  axisym- 
metric  nonlinear  finite  difference  calculations.  In  these  simulations,  spall  was 
observed  and  the  equivalent  seismic  sources  from  the  calculations  were 
modeled  as  a  point  explosion  plus  the  opening  and  closing  of  a  tension  crack 
(spall).  Such  a  model  is  illustrated  in  Figure  3.1.  In  this  section,  we  summarize 
results  of  calculations  that  attempt  to  use  these  models  to  predict  regional 
phases  such  as  Pg,  Lg.  and  Rg.  The  physically  reasonable  spall  models  are 
used  to  derive  consistent  models  for  the  equivalent  moment  tensor  representa¬ 
tions  for  the  seismic  source.  These  moment  tensor  sources  are  then  used  in 
conjunction  with  wavenumber  integration  Green’s  functions  for  the  response  of 
a  layered  earth  structure  to  compute  synthetic  regional  seismograms. 

As  part  of  this  work  we  present  a  new  derivation  for  the  approximate 
equivalence  of  the  surface  point  force  and  the  moment  tensor  representations 
of  a  buried  tension  crack.  The  approximation  is  tested  numerically  and  found  to 
be  adequate  for  tension  cracks  buried  less  than  1/4  of  the  shear  wave 
wavelength  for  Pg  and  Lg.  Proper  modeling  of  Rg  requires  wavelengths  longer 
than  eight  times  the  burial  depth.  In  addition,  a  smoothing  operator,  E(t),  is 
presented  that  results  from  a  model  of  finite  extent  of  the  spall  distributed  over 
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Figure  3.1.  The  axisymmetric  spall  above  the  explosion  is  modeled  as  a 
circular  horizontal  tension  crack  that  opens  and  closes  in  the 
vertical  direction.  The  spall  is  parameterized  by  the  source  depth 
hg,  the  crack  radius  a,  and  the  displacement  time  function  u(t,x). 
The  radiation  from  a  horizontal  tension  crack  is  equivalent  at  low 
frequencies  to  the  radiation  from  a  vertical  point  force,  F^,  at  the 
surface  proportional  to  the  second  time  derivative  of  8u.  The 
moment  tensor  source  representation  for  the  tension  crack  is 
proportional  to  5u.  A  is  the  area  of  the  crack,  p  is  the  density, 
and  X  and  p  are  the  Lame  parametes  for  the  medium. 
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the  tension  crack.  In  this  model,  motion  is  zero  at  the  crack  tip  and  spall 
motion  is  concave  upward.  Such  a  model  implies  that  the  spall  source  falls-oft 
at  high  frequencies  at  least  as  fast  as  Since  the  spall  source  is  propor¬ 
tional  to  a?  for  low  frequencies  and  falls-off  faster  than  at  high  frequencies, 
the  spall  contribution  will  always  be  a  narrowband  signal  compared  to  the  explo¬ 
sion  source  which  is  flat  at  low  frequencies  and  falls-off  as 

Figures  3.2  and  3.3  show  the  equivalent  force  time  histories  and  their 
spectra  for  teleseismic  P-wave,  Pn,  Pg,  and  Lg  slownesses.  Note  that  as  a 
consequence  of  the  finite  extent  of  spall  across  the  tension  crack,  the 
equivalent  seismic  source  depends  on  the  apparent  slowness  of  the  seismic 
phase.  Phases  with  slower  apparent  velocities  show  a  lower  corner  frequency 
and  appear  more  narrowband.  The  spall  momentum,  depth,  extent,  and  veloci¬ 
ties  are  given  by  the  Barker  and  Day  (1990)  parameterization  for  spall  observed 
in  the  finite  difference  simulations  for  a  125  KT,  680  meter  deep  explosion  in  a 
Shagan  River  model. 

Figure  3.4  shows  the  vertical  displacement  Lg  and  Pg  waveforms  predicted 
for  a  distance  of  300  km  from  the  explosion  and  the  spall  sources  in  Figures  3.2 
and  3.3.  The  seismic  sources  have  been  convolved  with  the  Green’s  functions 
for  an  explosive  source  and  those  for  a  tension  crack  to  predict  the  regional 
ground  motion  from  the  two  source  components.  Note  that  the  spall  source  has 
an  Lg  10  times  larger  than  the  explosion  point  source  Lg.  In  contrast  the  spall 
source  modulates  the  Pg  signal  interfering  with  the  pP  from  the  explosion  to 
complicate  the  Pg  waveform.  Spectra  are  shown  in  Figure  3.5  to  show  that  the 
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Figure  3.2.  The  equivalent  spall  source  (for  source  depth  of  680  m)  for  four 
different  apparent  phase  velocities  corresponding  to  teleseismic 
P,  Pn,  Pg,  and  Lg.  The  smoothing  operator  E{t)  simulates  the 
disk-like  nature  of  the  source  as  seen  in  the  far  field  by  a  wave  of 
a  given  slowness.  Slower  phases  are  more  sensitive  to  the  size 
of  the  spall  disk.  Because  the  spall  velocity  is  distributed  across 
the  tension  crack,  the  crack  opens  faster  than  it  closes  and  the 
seismic  source  is  asymmetrical.  The  slap  down  is  distributed 
over  a  longer  time  period  than  the  initial  spall  opening  and 
therefore  although  the  area  under  the  cunre  is  zero,  the  closing 
phase  has  lower  frequency  character  than  the  opening  phase. 
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EFFECTIVE  SPALL  SOURCE  FUNCTICV  SPEr[R 


FREQUENCY  (HZ) 


Figure  3.3.  Spectra  of  the  time  series  of  Figure  3.2.  The  finite  extent  of  the 
tension  crack  ensures  that  the  slower  phases  like  Lg  have  an 
equivalent  spall  source  that  falls*off  more  rapidly  than  the  higher 
phase  velocity  phases  such  as  Pg,  Pn  or  the  teleseismio  P*wave 
equivalent  spall  source. 
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Lg  WAVEFORMS  (H  =  680  M,  A  =  300  KM) 


SECONDS 


Pg  WAVEFORMS  (H=680  M.  A  =  300  KM) 


Figure  3.4.  Synthetic  Lg  and  Pg  waveforms  at  a  distance  of  300  km  for  the 
680  meter  depth  of  burial  125  KT  Shagan  River  explosion 
simulation.  The  spall  model  is  from  Barker  and  Day  (1990) 
based  on  parameterized  models  derived  from  a  nonlinear 
axisymmetric  finite  difference  calculation.  Note  that  the  Lg  spall 
signal  dominates  the  Lg  explosion  signal  while  the  Pg 
explosion+spall  signal  is  only  somewhat  larger  than  the  pure 
explosion  Pg  signal.  The  m|3(Lg)  from  the  spall  source  is  about 
6.2.  The  mjj(Lg)  from  the  pure  explosion  source  is  about  4.9. 
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Figure  3.5.  Spectra  of  the  time  series  in  Figure  3.4.  The  Lg  spall  signal  is 
about  20  times  larger  than  the  Lg  explosion  signal  at  1  Hz  while 
the  Pg  spall  signal  is  about  two  times  larger  than  the  pure 
explosion  Pg  signal  at  1  Hz.  The  Lg  spall  Signal  is  narrowband 
and  falls-off  rapidly  below  the  Lg  explosion  signal  above  3  Hz. 
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Lg  is  dominated  by  the  spall  source  for  frequencies  up  to  3  Hz.  The  spall 
source  provides  low  frequencies  below  1  Hz  and  the  spall+explosion  signal  is 
closer  to  the  average  of  the  two  sources  at  frequencies  above  1  Hz. 

We  conclude  that  the  parameterized  spall  models  of  Barker  and  Day 
(1990)  predict  significant  Lg  excitation  relative  to  the  point  explosion  part  for 
Eastern  Kazakhstan  crustal  models.  As  seen  in  Section  IV,  this  is  in  contrast  to 
crustal  structures  with  low  near  surface  velocities  and  high  near  surface 
attenuation  as  typical  of  basin  and  range  NTS  models.  Crustal  models  with 
high  near  surface  P-wave  velocities  exhibit  very  weak  Lg  excitation  from  point 
explosions.  The  shallow  tension  crack  source  in  these  models  Is  a  significant 
Lg  source  without  overly  dominating  the  Pg  waveform.  The  parameterized 
models  of  Barker  and  Day  (1990)  based  on  the  nonlinear  finite  difference  simu¬ 
lations  predict  an  (Lg)  =  6.2  for  the  125  KT  explosion.  In  contrast  the  point 
explosion  model  in  such  a  crustal  structure  predicts  an  (Lg)  =  4.9. 

While  there  are  other  hypotheses  for  the  excitation  of  Lg  from  explosions,  it 
is  clear  that  for  Eastern  Kazakhstan  structures,  the  spall  model  could  be  a 
significant  contributor  to  the  Lg  signal.  The  point  explosion  is  clearly  a  poor 
model  for  the  Lg  excited  by  an  explosion  in  these  structures.  Even  a  small 
amount  of  spall  will  overcome  the  Lg  excited  by  the  pure  explosion  part  of  the 
source. 
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Section  IV 

Analysis  of  Pg  and  Lg  Excitation 
by  Axisymmetric  Sources  in 
Layered  Anelastic  Crustal  Models 

Background 

Section  III  of  this  report  summarizes  work  from  McLaughlin,  et  al.  (1990) 
that  applies  the  Barker  and  Day  (1990)  parametric  seismic  source  models  sum¬ 
marized  in  Section  II  and  shows  they  imply  spall  should  be  a  significant  source 
of  Lg  in  addition  to  the  pure  explosive  source.  In  order  to  put  this  conclusion 
into  perspective,  we  report  in  this  section  on  numerical  calculations  that  demon¬ 
strate  differences  between  Pg  and  Lg  excitation  by  point  explosion,  spall,  and 
CLVD  sources  in  layered  earth  structures  representative  of  the  basin  and  range 
and  the  Eastern  Kazakhstan  regions.  The  crustal  structures  for  these  two 
regions  are  very  different  and  lead  to  significant  differences  in  theoretical  Pg 
and  Lg  excitation  by  simple  point  sources.  Generalizations  about  relative  exci¬ 
tation  of  regional  phases  by  different  source  types  in  one  structure  often  do  not 
apply  in  the  other  structure. 

In  addition  to  the  simple  point  explosion  source  we  examine  the  Pg  and  Lg 
excitation  from  axisymmetric  deviatoric  sources  suggested  by  the  spall  models 
of  Barker  and  Day.  There  have  been  previous  studies  that  compare  the  excita¬ 
tion  of  Pg  and  Lg  by  shallow  point  explosions  for  purposes  of  examining  the 
theoretical  basis  of  explosion  and  earthquake  discrimination,  see  for  example, 
Campillo,  et  al.  (1984),  Bennett,  et  al.  (1987),  or  Lilwall  (1988).  This  short 
report  however,  differs  from  previous  work  in  that  it  examines  some  of  the 
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nonspherically  symmetric  components  of  the  contained  underground  explosion. 
Furthermore,  the  use  of  a  wavenumber  integration  algorithm  permits  us  to 
model  the  attenuation  of  the  crust  in  a  more  exact  manner  than  the  modal  sum¬ 
mation  or  discrete  wavenumber  methods  used  by  some  earlier  workers. 

Any  nonlinear  free-surface  interaction  must  be  accompanied  by  shallow 
deviatoric  sources.  Therefore,  the  relative  importance  of  these  secondary  non¬ 
linear  sources  for  the  excitation  of  regional  phases  is  important  for  stable  yield 
estimators  as  well  as  for  discrimination  between  earthquakes  and  explosions. 
Spall  is  the  simplest  nonlinear  free-surface  interaction  that  is  commonly 
observed,  but  by  no  means  are  the  shallow  nonlinear  interactions  limited  to  just 
spall. 

When  we  consider  the  factors  that  may  cause  the  explosion  source  to  devi¬ 
ate  from  spherical  symmetry,  it  is  clear  that  the  presence  of  the  free-surface, 
the  vertical  gradient  in  confining  stress,  and  the  vertical  variation  of  material 
properties  are  the  most  significant.  All  of  these  effects  may  be  considered  in  a 
two-dimensional  axisymmetric  geometry  (cylindrical  coordinates)  where  material 
properties  are  a  function  of  only  depth  (z). 

For  simplicity  we  restrict  ourselves  to  the  study  of  axisymmetric  sources 
although  there  is  ample  evidence  that  explosions  at  the  Nevada  Test  Site  (NTS) 
radiate  long  period  Love  waves  (Toksoz  and  Kehrer,  1972),  long  period  SH 
body  waves  (Wallace,  et  al.,  1985),  as  well  as  short  period  SH  body  waves 
(McLaughlin,  et  al.,  1983;  Stump  and  Johnson,  1984;  Johnson,  1988).  Explo¬ 
sions  at  the  Eastern  Kazakhstan  Test  Site  (EKTS)  are  known  to  exhibit  tectonic 
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release  in  long-period  surface  waves  (Day  and  Stevens,  1986)  and  recently 
anomalous  SH  phases  have  been  observed  at  regional  distances  (Priestley, 
et  al. ,  1990).  Although  the  origins  of  SH  waves  from  explosions  are  still  under 
investigation,  no  researcher  has  found  a  consistant  correlation  between  long- 
period  tectonic  release  and  short-period  seismic  signals.  Johnson  (1988)  found 
that  the  non-diagonal  parts  (non-axisymmetric  components)  of  the  moment  ten¬ 
sor  (0.25-5  Hz)  were  small  and  the  largest  deviation  from  sphericity  was  the 
variation  in  versus  Hy  appears  that  the  most  important  deter¬ 

ministic  nonspherical  aspects  of  the  explosion  source  can  be  characterized  by 
axisymmetric  sources. 

We  consider  three  generic  axisymmetric  moment  tensor  sources 
Mj,y=/lf„=/Wy2=0):  the  explosion,  the  compensated  linear  vector  dipole  (CLVD), 
and  the  horizontally  oriented  tension  crack.  The  linear  combination  of  any  two 
of  these  sources  will  serve  to  represent  any  axisymmetric  first  order  moment 
tensor  source.  The  explosion  moment  tensor, 
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where  A/f/  is  the  isotropic  explosion  moment,  Mqlvd  is  the  CLVD  moment,  M^c 
is  the  tension  crack  moment,  and  5i7  is  the  displacement  separation  across  the 
tension  crack  with  area  A.  The  tension  crack  has  been  used  as  a  model  for 
spall  (Day,  et  al,  1983),  although  general  axisymmetric  nonlinear  free-surface 
interactions  may  result  in  more  general  linear  combinations  of  any  of  these 
three  generic  sources.  It  is  unlikely,  however,  that  any  nonlinear  deformation 
would  occur  that  did  not  exhibit  some  deviatoric  component  since  the  normative 
failure  mode  in  rocks  is  shear  failure. 

For  the  purposes  of  computing  far-field  regional  seismograms  we  se  a 
wavenumber  integration  algorithm  (Apsel  and  Luco,  1983)  that  computes  the 
complete  P-SV  wavefield  in  an  attenuating  layered  earth  model.  This  method 
has  been  found  to  be  superior  to  modal  summation  for  regional  seismogram 
synthesis  when  attenuation  is  a  significant  aspect  of  the  propagation  (Day, 
et  al.,  1989). 

Several  crustal  models  are  considered  to  give  a  flavor  for  the  potential  vari¬ 
ation  of  regional  phase  excitation  from  model  to  model  but  the  models  are  not 
considered  to  be  exhaustive.  Tabulations  of  the  models  are  found  at  the  end  of 
this  section  in  Tables  4.1  through  4.8.  Two  basic  sets  of  models  are  con¬ 
sidered.  First,  models  with  a  high  free-surface  velocity  (P-wave  velocity  >  4000 
m/s)  are  considered.  Next  some  models  with  low  free-surface  velocities  (P- 
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wave  velocity  <  4000  m/s)  and  hence  high  velocity  gradients  in  the  upper  crust 
are  considered.  We  are  concerned  with  evaluating  the  relative  excitation  of  Pg 
and  Lg  as  a  function  of  frequency  in  the  short-period  band  (0.25  to  5  Hz)  for 
shallow  axisymmetric  moment  tensor  sources.  An  important  aspect  of  many  of 
these  models  is  that  attenuation  (1/Q)  is  higher  in  the  upper  few  kilometers  of 
the  crust.  This  aspect  of  a  model  serves  to  attenuate  high  frequency  Rg  and 
other  waves  that  are  confined  to  the  near-surface  layers.  This  is  particularly 
important  in  models  with  low  near-surface  velocities. 

High  Near-Surface  Velocity  Models 

Lilwall  (1988)  presented  a  theoretical  analysis  of  the  discrimination 
between  earthquakes  and  explosions  using  explosion  and  earthquake  sources 
located  in  a  crustal  structure  with  an  upper  1000  m  thick  layer  of  4000  m/s  P- 
wave  velocity.  He  points  out  that  in  such  a  high  P-wave  velocity  structure  that 
Lg  amplitude  is  largely  due  to  a  ‘nongeometric’  S  wave  referred  to  as  S*.  The 
curved  compressional  wavefront  from  the  point  explosion  source  produces  a 
frequency  dependent  SV  reflection  that  is  seen  at  regional  distances  as  Lg.  We 
examine  Lilwall’s  suggestion  that  the  S*  phase  is  largely  responsible  for  the  Lg 
from  shallow  explosions.  If  true,  then  the  low-frequency  character  of  the  S* 
phase  would  provide  a  theoretical  basis  for  the  discrimination  between  explo¬ 
sions  and  earthquakes. 

Figure  4.1  shows  a  suite  of  Green’s  functions  for  explosion  sources  at 
depths  of  100,  300,  500,  700,  900,  1100  and  1300  meters  for  a  surface  vertical 
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Figure  4.1 .  Explosion  Green’s  functions  (0-5  Hz)  at  a  distance  of  300  km  for 
source  depths  of  100  to  1300  m.  Source  is  a  1  Nt-m  step 
function  and  vertical  velocity  is  in  units  of  10"^^  m/s.  Pg,  Lg  and 
Rg  wavepackets  are  labeled.  Note  the  decay  of  the  excitation  of 
Lg  with  increasing  depth  of  the  source. 
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receiver  at  a  range  of  300  km.  The  structure  is  derived  from  Stevens  (1986) 
and  has  a  surface  P-wave  velocity  of  5020  m/s.  The  Green’s  functions  are 
broadband  (0-5  Hz)  velocity  responses  (units  of  m/s)  to  a  step  function  explo¬ 
sion  source  (units  of  Nt-m).  The  Pg  amplitude  is  roughly  constant  modulated  by 
variation  in  the  P+pP  interference  with  depth.  The  broadband  Lg  amplitude 
decreases  with  depth  as  expected  if  the  primary  source  for  SV  energy  were  the 
conversion  of  a  curved  P  wavefront  at  the  free-surface.  The  Lg  is  also  fre¬ 
quency  dependent.  The  spectral  ratios  of  Lg/Pg  are  shown  in  Figure  4.2  as  a 
function  of  frequency.  At  0.25  Hz,  the  Lg  is  larger  than  the  Pg  for  all  depths 
and  larger  for  the  100  meter  source  than  the  1300  meter  deep  source.  The 
variance  between  the  100  and  1300  meter  deep  sources  is  never  greater  than  a 
factor  of  5  and  is  much  less  clear  at  frequencies  above  2  Hz  than  below  2  Hz. 
Note  that  at  1  Hz,  all  explosion  sources  exhibit  Lg/Pg  ratios  less  than  1 . 

At  the  same  scale  as  Figure  4.1,  a  suite  of  tension  crack  Green’s  functions 
for  the  same  source  depths  are  shown  in  Figure  4.3.  For  the  same  moment, 
the  tension  crack  Pg  amplitudes  are  about  1/2  the  explosion  Pg  amplitudes 
while  the  tension  crack  Lg  amplitudes  are  several  times  larger  than  the  explo¬ 
sion  Lg  amplitudes.  The  deviatoric  part  of  the  tension  crack  source  provides  a 
direct  source  of  SV  energy  to  propagate  at  Lg  phase  velocities.  An  investiga¬ 
tion  of  the  frequency  content  of  the  Lg/Pg  ratio  shows  that  the  tension  crack 
produces  more  Lg  than  an  explosion  source  over  a  broad  bandwidth.  The 
discrepancy  is  so  large  that  small  amounts  of  a  tension  crack  source  will  pro¬ 
duce  Lg  comparable  to  the  pure  explosion. 
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=  300  km  in  a  high  velocity  stmcture  (Stevens,  1986). 
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Figure  4.3.  Tension  crack  Green’s  functions  (0-5  Hz)  at  a  distance  of  300  km 
for  source  depths  of  100  to  1300  m.  Source  is  a  1  Nt-m  step 
function  and  vertical  velocity  is  in  units  of  10*^®  m/s.  Pg,  Lg,  and 
Rg  wavepackets  are  labeled.  Note  that  in  comparison  to 
Figure  4.1 ,  the  Lg  is  larger  and  does  not  decay  with  increasing 
depth  of  the  source. 
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Figures  4.4  and  4.5  show  Green’s  functions  at  a  receiver  distance  of  256 
km  for  an  explosion  and  a  CLVD  source  at  a  depth  of  500  meters  as  computed 
for  a  suite  of  models  proposed  for  the  region  around  the  EKTS  (see  Tables 
4. 1-4.7).  The  models  have  varying  crustal  thickness,  Pn  velocities,  attenuation, 
and  crustai  gradients.  Three  of  the  models  have  crustal  low-velocity  zones.  All 
models  have  surface  P-wave  velocities  greater  than  5000  m/s.  All  of  the 
models  exhibit  very  low  Lg/Pg  levels  from  the  explosion  and  larger  Lg/Pg  ratios 
for  the  CLVD  source.  It  is  clear  that  these  results  are  not  sensitive  to  the 
details  of  the  crustal  models. 

Another  way  to  consider  this  result  is  that  a  1 0%  asymmetry  in  the  moment 
tensor  source  can  result  in  as  much  Lg  as  the  spherical  source  would  produce. 
Any  nonspherical  expansion  or  spall  will  result  in  a  CLVD  component  superim¬ 
posed  on  an  explosion  source.  The  Lg  excited  by  the  deviatoric  CLVD  source 
will  mask  the  explosion  produced  S*. 

Low  Near-Surface  Velocity  Models 

Crustal  models  with  a  low  near-surface  velocity  will  invariably  have  a 
higher  velocity  gradient  in  the  upper  crust.  This  low  velocity  channel  traps 
waves  with  low  phase  velocity  excited  by  a  shallow  source.  These  low  velocity 
near-surface  layers  are  not  contiguous  at  NTS.  However,  explosions  at  NTS 
initially  excite  waves  in  these  low  velocity  layers  and  the  near-source  structure 
determines  the  character  of  the  energy  radiated  into  the  crustal  waveguide.  To 
properly  treat  where  this  shallow  propagating  energy  goes  when  these  waves 
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Figure  4.4.  Explosion  Green’s  functions  (0-4  Hz)  at  a  distance  of  256  km  for 
a  suite  of  different  crustal  structures.  The  source  is  at  a  fixed 
depth  of  500  m.  Note  the  weak  Lg  excitation  for  all  models.  The 
top  model  (labeled  Stevens)  has  a  low  Q  layer  near  the  surface 
that  attenuates  the  high-frequency  Rg. 
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Figure  4.5.  CLVD  Green’s  functions  {0-4  Hz)  at  a  distance  of  256  km  for  a 
suite  of  different  crustal  structures.  The  source  is  at  a  fixed 
depth  of  500  m.  Note  the  strong  Lg  excitation  for  all  models. 
The  top  model  (labeled  Stevens)  has  a  low  Q  layer  near  the 
surface  that  attenuates  the  high-frequency  Rg.  Traces  have 
been  clipped  at  +/-0.0015  in  order  to  display  them  at  the  same 
scale  as  Figure  4.3. 


51 


SSS-TR-90-11536 


impinge  on  edges  of  these  low  velocity  channels  is  a  scattering  problem  beyond 
this  current  work.  McLaughlin,  et  al.  (1987)  suggest  that  some  of  this  energy  is 
ultimately  scattered  into  the  crustal  waveguide  as  Lg,  so  the  analysis  presented 
here  may  be  modified  by  the  scattering  hypothesis  for  the  generation  of  Lg.  For 
the  purposes  of  this  report  we  simply  wish  to  show  that  the  relative  Pg  and  Lg 
excitation  by  explosive  and  CLVD  sources  in  a  basin  and  range  structure  with 
low  velocity  and  low  Q  near-surface  layers  is  quite  different  from  the  excitation 
in  the  high  near-surface  velocity  models  previously  discussed. 

BR3  (derived  from  BR2,  McLaughlin,  et  al.,  1988)  is  a  layered  crustal  and 
upper  mantle  model  with  constant  Q,  designed  to  fit  the  propagation  of  Pn,  Pg 
and  Lg  waveforms  in  the  1  Hz  bandwidth  for  explosions  at  NTS  and  stations  in 
the  basin  and  range.  The  model  correctly  predicts  the  spatial  attenuation  of 
these  phases  in  the  regional  distance  range.  At  frequencies  substantially  above 
1  Hz,  the  model  progressively  fails  to  model  the  Lg  attenuation  properly 
because  apparent  Lg  Q(f}  increases  with  increasing  frequency  and  BR3  is  a 
constant  Q  model.  BR3  is  characterized  by  low  velocity  and  low  Q  layers  at  the 
surface  with  velocities  and  Q  increasing  with  depth.  The  model  was  chosen  not 
to  optimize  the  fit  for  any  one  path  but  to  represent  an  average  over  paths. 

Green’s  functions  (0-4  Hz)  for  a  distance  of  320  km  are  shown  in  Figure 
4.6  for  explosion  and  CLVD  sources  at  depths  of  250  and  550.  In  contrast  to  a 
similar  comparison  that  is  made  in  Figure  4.4  and  4.5  for  the  high  velocity  struc¬ 
ture,  we  see  that  in  the  BR3  structure,  the  CLVD  does  not  excite  much  more  Lg 
than  the  explosion.  The  Eastern  Kazakhstan  and  BR3  models  lead  to  very 
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Figure  4.6.  Explosion  and  CLVD  Green’s  functions  (0-4  Hz)  at  a  distance  of 
320  km  in  the  BR3  model.  Sources  are  either  250  or  550  m 
deep.  Note  that  the  peak  Lg  amplitude  from  the  250  m  CLVD  is 
about  the  same  as  the  Lg  amplitude  from  the  550  meter 
explosion.  The  CLVD  sources  excite  higher-modes  that  arrive 
superimposed  over  the  LR. 
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different  conclusions  regarding  the  Lg  excitation  from  CLVD  and  explosion 
sources. 

In  the  BR3  model,  Lg  from  the  CLVD  is  more  complicated  than  the  Lg  from 
the  explosion  sources.  It  appears  that  the  CLVD  excites  higher  modes  that 
interfere  with  the  early  portions  of  the  Rayleigh  wave.  However,  the  peak  Lg 
amplitude  from  the  CLVD  source  is  not  significantly  larger  than  the  Lg  ampli¬ 
tudes  due  to  the  explosion  sources.  In  a  stmcture  such  as  BR3,  the  CLVD  and 
explosion  sources  yield  similar  Lg/Pg  ratios.  Therefore,  in  the  BR3  structure, 
the  Lg/Pg  ratio  is  invariant  with  respect  to  the  three  generic  axisymmetric 
sources. 

Conclusions 

The  two  sets  of  structures  considered,  show  substantially  different  relative 
responses  to  shallow  explosion  and  CLVD  source  components.  Explosions  in 
structures  with  near-surface  high  velocities  couple  little  SV  energy  into  the  Lg 
waveguide  while  a  shallow  CLVD  component  produces  much  larger  SV  energy 
which  couples  into  the  Lg  waveguide.  Consequently  the  CLVD  component  of 
the  Lg  in  the  high  velocity  structure  may  easily  dominate  the  explosion  com¬ 
ponent  of  the  Lg.  In  stark  contrast,  the  shallow  explosion  and  CLVD  sources  in 
the  low  velocity  structure  produce  comparable  Pg/Lg  ratios.  It  appears  that  the 
high  attenuation  of  the  near-surface  low-Q  layer  attenuates  shear  waves 
trapped  in  the  near-surface  layers  while  the  SV  energy  from  the  explosion  and 
CLVD  components  available  to  the  Lg  waveguide  is  more  nearly  equal  in  the 
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low-velocity  structure.  However,  the  CLVD  Lg  signal  in  the  low-velocity  struc¬ 
ture  is  more  complex  presumably  because  it  excites  higher  modes  in  the  sub  1 
Hz  frequency  band.  Patton  (1988)  has  used  these  higher  modes  to  model  the 
CLVD  components  of  NTS  explosions.  Because  the  explosion  and  CLVD 
sources  produce  roughly  equal  Lg/Pg  ratios  a  procedure  similar  to  Patton’s  is 
required  to  decipher  the  more  complex  nature  of  the  source.  The  indeterminacy 
of  Lg-to-Pg  ratios  at  NTS  is  demonstrated  in  Section  V  of  this  report  using  two 
co-located  NTS  events  recorded  at  a  common  station. 

It  is  likely  that  any  underground  explosion  contains  some  asymmetries. 
The  simplest  first  order  moment  tensor  source  which  models  asymmetry  in  the 
explosion  source  is  a  CLVD.  Given  the  sensitivity  of  the  Lg  excitation  in  a 
high-velocity  crustal  structure  to  the  CLVD  component  it  seems  unlikely  that  the 
observed  regional  Lg  is  generated  by  an  explosion  S*  mechanism.  Whether 
this  asymmetry  is  due  to  spall  or  the  vertical  gradient  in  the  material  properties 
and  confining  pressure  is  immaterial  to  the  argument.  Since  the  Lg  excitation 
by  a  CLVD  source  is  about  an  order  of  magnitude  greater  than  by  an  explosion 
source,  only  a  10%  CLVD  component  is  required  to  generate  as  much  Lg  as 
the  pure  explosion. 

Johnson  (1988)  has  found  from  the  inversion  of  near-field  accelerograms  at 
NTS  that  the  component  of  the  moment  tensor  was  out  of  phase  with  the 
Mxx  and  Myy  components  introducing  a  CLVD  source  that  was  consistent  with 
spall.  Patton  (1988)  examining  the  same  event,  HARZER,  found  similar  evi¬ 
dence  for  the  CLVD  component  of  the  source  while  modeling  higher-mode 
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components  of  the  Lg  wavetrain.  In  McLaughlin,  et  at.  (1990)  (summarized  in 
Section  III  of  this  report)  we  show  that  the  parametric  spall  models  of  Barker 
and  Day  (1990)  (summarized  in  Section  II  of  this  report)  predict  (Lg)  values 
within  reasonable  bounds  of  observation.  Whether  the  axisymmetric  asym¬ 
metries  in  the  equivalent  seismic  source  are  due  to  spall  or  other  nonlinear 
processes,  it  is  clear  that  these  nonlinear  processes  will  dominate  the  excitation 
of  Lg  in  crustal  models  with  near-surface  high  velocities. 

Given  that  Lg  is  known  to  be  a  very  stable  yield  estimator,  the  question 
arises  as  to  whether  asymmetries  in  the  explosion  source  could  be  responsible 
for  the  reproducible  excitation  of  Lg  energy  proportional  to  the  yield.  Perhaps, 
for  explosions  detonated  in  a  fairly  uniform  brittle  material  at  nearty  constant 
scaled  depth,  the  spall  momentum  is  a  good  estimate  of  yield.  Other  mechan¬ 
isms  also  may  be  at  work  to  stabilize  the  process  of  Lg  excitation.  Scattering  of 
Rg  (and  shallow  P-SV  modes)  may  contribute  to  the  regional  Lg  signal.  Shal¬ 
low  explosions  are  efficient  Rg  sources.  Rg  may  be  attenuated  by  scattering  as 
well  as  by  intrinsic  attenuation  in  the  upper  crust  (McLaughlin  and  Jih,  1987;  Jih 
and  McLaughlin,  1988).  Since  scattering  occurs  over  a  large  volume  of  crust 
there  would  be  a  statistical  averageing  that  may  stabilize  the  Lg  signal.  Untill  we 
have  a  better  understanding  of  these  processes,  our  understanding  of  the 
Lg:yield  relation  will  be  unsatisfactory. 
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Tabulated  Crustal  Models 

The  model  in  Table  4.1  is  derived  from  a  model  Stevens  (1986)  based  on 
surface  wave  studies  with  some  modification  due  to  DSS  reports  for  the  region. 
The  Q’s  have  been  increased  in  the  crust  so  the  average  crustal  Q  is  represen¬ 
tative  of  reported  Q’s  for  the  region.  The  model  of  Leith  (1987)  was  compiled 
from  Soviet  literature.  The  Priestley  et  al.  (1 988)  models  are  based  on  telese- 
ismic  P-wave  receiver  functions  for  the  stations  BAY  and  KKL.  The  models  from 
Thurber,  et  al.  (1989)  were  used  to  locate  quary  blasts  in  the  region  and  are 
derived  to  some  extent  from  the  Priestley,  et  al.  (1988)  models.  For  most 
models  (3  =  c(/'/3,  and  p  =  770  +  0.32  a.  Where  was  not  specified,  a  uni¬ 
form  value  of  367  was  used  as  reported  by  Priestley,  et  al.  (1990).  Therefore 
these  tabulations  represent  models  actually  used  in  calculations  and  not  neces¬ 
sarily  the  exact  models  advocated  by  the  cited  authors.  Model  BR3  is  a 
modified  version  of  BR2  from  McLaughlin,  et  al.  (1988)  with  a  lower  near¬ 
surface  Q  compensated  by  a  higher  Q  in  the  lower  crust. 


57 


sss-TR-go-nsae 


Table  4.1. 

Derived  from  Stevens  (1986) 

1 

2000 

5020 

2790 

2700 

100 

OO 

1000 

5400 

3000 

2700 

150 

oo 

2488 

5900 

3300 

2700 

200 

oo 

10976 

6100 

3400 

2700 

600 

oo 

5488 

6308 

3541 

2702 

525 

oo 

5488 

6597 

3703 

2807 

500 

oo 

5488 

6736 

3781 

2858 

450 

oo 

5564 

6782 

3807 

2875 

400 

oo 

6504 

6795 

3814 

2879 

350 

oo 

8006 

8147 

4573 

3372 

179 

oo 

9359 

8138 

4568 

3369 

167 

oo 

10940 

8106 

4550 

3358 

159 

oo 

12780 

8065 

4527 

3343 

153 

oo 

14950 

8047 

4517 

3336 

150 

oo 

17470 

8070 

4530 

3345 

148 

oo 

20420 

8117 

4556 

3361 

148 

oo 

23880 

8154 

4577 

3375 

147 

oo 

27910 

8161 

4581 

3378 

147 

oo 

32630 

8145 

4572 

1372 

146 

oo 

38140 

8120 

4558 

3363 

146 

oo 

OO 

8101 

4547 

3356 

145 

oo 

Table  4.2.  Leith  (1987)  DSS  Model 
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Table  4.3.  Priestley,  et  al. 

(1988)  BAY  Model  J 

h(m) 

a(m/s) 

1000 

5320 

3071 

2472 

367 

OO 

3000 

5960 

3441 

2677 

367 

oo 

6000 

5660 

3267 

2581 

367 

oo 

5000 

6010 

3469 

2693 

367 

oo 

10000 

6290 

3631 

2782 

367 

oo 

5000 

6650 

3839 

2898 

367 

oo 

15000 

6890 

3977 

2974 

367 

oo 

oo 

8330 

4809 

3435 

367 

oo 

Table  4.4.  Preistlev.  e  al. 

(1988)  KKL  ModelJ  1 

h(m) 

a(m/s) 

p(kg/m^) 

_Qu 

Ok 

1000 

5190 

2996 

2430 

367 

OO 

4000 

5850 

3377 

2642 

367 

oo 

2500 

5770 

3331 

2616 

367 

OO 

2500 

5570 

3215 

2552 

367 

oo 

2500 

5900 

3406 

2658 

367 

oo 

5000 

6190 

3573 

2750 

367 

oo 

10000 

6520 

3764 

2856 

367 

oo 

15000 

7040 

4064 

3022 

367 

oo 

10000 

7470 

4312 

3160 

367 

oo 

5000 

7800 

4503 

3266 

367 

oo 

OO 

8100 

4676 

3362 

367 

oo 

Table  4.5.  Priestley,  et  al.  (1988)  KKL  Model  S 


h(m) 

a(m/s) 

P(m/s) 

WBSSSMKSMEM 

1000 

5140 

2967 

2414 

367 

OO 

4000 

5910 

3412 

2661 

367 

oo 

7500 

5510 

3181 

2533 

367 

oo 

2500 

6030 

3481 

2699 

367 

oo 

12500 

6510 

3758 

2853 

367 

oo 

12500 

7080 

4087 

3035 

367 

oo 

2500 

7290 

4208 

3102 

367 

oo 

15000 

7490 

4324 

3166 

367 

oo 

OO 

8100 

4676 

3362 

367 

oo 
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Table  4.6.  Thurber,  Given,  and  Berger  (1989)  Model  A 


h(m) 

a{m/s) 

3(m/s) 

D(ka/m^) 

Ok 

5000 

5400 

3050 

2498 

367 

OO 

5000 

6150 

3500 

2738 

367 

oo 

10000 

3600 

2802 

367 

oo 

10000 

6550 

3700 

2866 

367 

oo 

10000 

6750 

3850 

2930 

367 

oo 

10000 

6950 

3950 

2994 

367 

oo 

oo 

8200 

4650 

3394 

367 

oo 

Table  4.7. 

Thurber. 

Given,  and  Beraer  (19891  Model  B 

h(m) 

a(m/s) 

P(m/s) 

HSSSSKI 

mm 

5000 

5400 

3300 

2498 

367 

OO 

5000 

6150 

3400 

273o 

367 

oo 

10000 

6350 

3500 

2802 

367 

oo 

10000 

6550 

3700 

2866 

367 

oo 

10000 

6750 

4100 

2930 

367 

oo 

10000 

6950 

4300 

2994 

367 

oo 

OO 

8200 

4700 

3394 

367 

oo 
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Table  4.8.  BR3,  NTS  -  Basin  &  Ranqe  Model 

h(m) 

a(m/s) 

p(kg/m^) 

Qu 

Ok 

500 

3000 

1500 

2000 

50 

oo 

500 

3800 

2200 

2200 

75 

oo 

2000 

5800 

3400 

2800 

125 

oo 

3000 

6000 

3460 

2900 

125 

oo 

14000 

6300 

3640 

2900 

150 

oo 

4000 

6500 

3750 

2900 

150 

oo 

4000 

6800 

3930 

3000 

150 

oo 

6000 

6900 

3980 

3000 

150 

oo 

4000 

7800 

4500 

3200 

JOO 

oo 

4000 

7840 

4500 

3200 

200 

oo 

4000 

7880 

4500 

3200 

200 

oo 

4000 

7920 

4500 

3200 

200 

oo 

4000 

7960 

4500 

3200 

200 

oo 

2000 

8000 

4500 

3200 

200 

oo 

4000 

8050 

4500 

3200 

200 

oo 

10000 

7900 

4400 

3200 

75 

oo 

10000 

7800 

4200 

3200 

75 

oo 

20000 

7800 

4200 

3200 

75 

oo 

20000 

7900 

4500 

3200 

75 

oo 

60000 

8000 

4550 

3300 

75 

oo 

etc 

8100 

4600 

3400 

150 

oo 
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Section  V 

An  Analysis  of  Two  Co-Located  NTS  Events 
with  Different  Scaled  Depths  of  Burial 

It  has  been  suggested  that  nonlinear  free-surface  effects,  such  as  spall, 
may  be  efficient  radiators  of  SV  radiation  and  hence  regional  Lg.  Nonlinear 
free-surface  effects  are  expected  to  scale  with  depth  in  such  a  way  that  events 
at  greater  scaled  depth  should  exhibit  less  spall  and  hence  excite  less  SV  radi¬ 
ation  from  spall.  Testing  such  a  hypothesis  by  the  comparison  of  two  similar 
sized  events  at  significantly  different  scaled  depths  is  difficult  at  NTS.  This  is 
because  of  the  variability  of  the  geologic  structure  and  because  coupling  is 
known  to  vary  across  the  water  table  at  NTS.  Large  events  at  normal  scaled 
depth  of  burial  or  small  events  at  large  scaled  depth  of  burials  are  usually  below 
the  water  table  while  smaller  events  at  normal  depth  of  burial  are  above  the 
water  table.  Also,  unless  two  events  are  close  together  differences  in  path 
effects  can  make  detailed  comparisons  difficult. 

The  underground  explosions,  Buteo  and  Duryea,  were  two  events 
detonated  in  the  mid  I960’s  on  Pahute  Mesa  at  quite  different  scaled  depths 
but  virtually  the  same  location  (see  Table  5.1).  The  water  table  at  this  location 
was  determined  to  be  at  a  depth  of  660  m,  so  Duryea  was  placed  above  the 
saturated  zone  in  rhyolite  while  Buteo  was  placed  in  saturated  tuff.  Duryea  had 
an  announced  yield  of  65  KT  and  was  therefore  at  a  scaled  depth  of  135 
which  is  within  the  usual  range  of  scaled  depths  for  testing  at  NTS. 
Buteo  was  substantially  smaller  and  therefore  over  buried.  As  we  will  see 
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below,  based  on  the  1  Hz  amplitude  ratios  of  Pg  and  Lg  at  KN-UT,  Buteo  was 
between  20  and  30  times  smaller  than  Duryea  and  therefore  was  at  a  scaled 
depth  of  at  least  450  or  deeper.  Nuttli  (1986)  included  Duryea 

{mt,{Lg)  =  5.66)  with  other  well  coupled  events  at  NTS  and  found  that  the  Lg 
excitation  from  Duryea  was  consistent  with  other  events  in  water-saturated  rock 
at  NTS.  Therefore,  although  the  two  events  are  not  of  similar  size  and  above 
and  below  the  water  table,  the  coupling  appears  to  be  excellent  for  both  events 
and  they  are  nearly  co-located  with  a  large  difference  in  scaled  depth  of  burial. 


Table  5.1.  Event  Information 

Event 

Date 

O.T. 

Lat. 

HEBI3H 

Medium 

Buteo 

12  May  1965 

18:15:00.1 

37.24 

-116.43 

696. 

tuff 

Duryea 

14  April  1966 

14:13:43.1 

37.24 

-116.43 

544. 

rhyolite 

Since  the  two  events  are  virtually  at  the  same  location,  comparisons  of  Pg 
and  Lg  at  a  common  seismic  station  can  serve  to  test  the  hypothesis  that  varia¬ 
tions  in  scaled  depth  cause  variations  in  the  relative  Pg  to  Lg  excitation. 
Although  the  different  yields  imply  different  explosion  sources,  the  ratios  of 
Lg/Pg  from  two  point  explosion  sources  so  close  together  should  be  the  same 
recorded  at  the  same  regional  station.  Few  on  scale  recordings  of  these  two 
events  are  available  at  common  stations:  however,  the  LRSM  station  KN-UT 
(Kanab  Utah)  at  a  distance  of  321  km  did  record  these  two  events  with  good 
signal-to-noise  in  the  short  period  band. 
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KN-UT  seismograms  are  shown  for  the  two  events  in  Figure  5.1.  Because 
highpass  filtered  records  of  the  two  events  showed  that  the  Lg  above  5  Hz  was 
primarily  Pg  coda,  we  use  bandpass  filtered  estimates  of  the  signal  strength 
rather  than  spectral  estimates.  In  this  way,  the  Lg  excitation  is  only  measured  if 
energy  can  be  seen  above  the  coda  of  the  preceding  phase.  Bandpass  filters 
were  applied  to  the  vertical  and  radial  components  and  the  rectified  seismogram 
was  smoothed  with  a  2  second  moving  average.  The  peak  amplitude  of  the 
smoothed  envelope  was  taken  as  a  measure  of  the  excitation  of  Pg  and  Lg  in 
each  frequency  bandpass.  The  procedure  was  repeated  for  radial  component 
of  motion.  Ratios  of  Lg/Pg  and  Lg/Lg  were  found  to  be  the  same  using  both  the 
vertical  and  radial  components  so  the  results  shown  are  based  on  the  average 
of  the  vertical  and  radial  components  of  motion. 

The  upper  portion  of  Figure  5.2  shows  the  ratio  of  (Duryea  Pg)/(Buteo  Pg) 
and  (Ouryea  Lg)/(Buteo  Lg).  Note  that  both  phases  show  a  systematic  decline 
in  the  ratio  from  0.25  Hz  to  5  Hz.  At  1  Hz,  Duryea  is  about  20  times  larger  than 
Buteo  while  at  5  Hz,  Duryea  is  only  about  5  to  6  times  larger  than  Buteo.  This 
is  consistent  with  Buteo's  higher  corner  frequency  since  it  was  a  smaller  yield 
and  located  at  a  greater  depth.  The  Lg  window  is  dominated  by  Pg  coda  above 
5  Hz  and  no  Lg/Lg  ratios  are  shown  for  frequencies  above  5  Hz.  The  (Duryea 
Lg)/(Duryea  Pg)  and  (Buteo  Lg)/(Buteo  Pg)  ratios  are  shown  in  the  lower  por¬ 
tion  of  Figure  5.2.  The  relative  excitation  of  Lg/Pg  is  quite  consistent  for  the 
two  events  up  to  5  Hz.  Below  1  Hz,  the  Lg/Pg  ratio  is  about  3  to  2  while  above 
1  Hz,  the  Lg/Pg  ratio  is  about  1  to  4  for  both  events.  Upper  bounds  to  the 
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DURYEA  AT  KN-UT  SPZ,  A-321km 


BUTEO  AT  KN-UT  SPZ,  A  =  321km 


Time  (Sec) 


Figure  5.1 .  Duryea  and  Buteo  recorded  on  the  short  period  vertical  channel 
at  LRSM  station  KN-UT  (Kanab,  Utah).  Pn,  Pg  and  Lg  phases 
are  indicated  on  the  figure.  The  similarity  of  the  two  records  is 
clear  with  the  Buteo  exhibiting  a  somewhat  higher  frequency 
signal  than  Duryea  due  to  source  scaling. 
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DURYEA/BUTEO  RATIO  AT  KN-UT 


0  1  2  3  4  5  6  7 

FREQUENCY  (Hz) 


LG/PG  RATIO  AT  KN-UT 


o  ^ 

t  ’  UPPER  BOUND 


FREQUENCY  (Hz) 

Figure  5.2.  Ratios  of  bandpass  filtered  KN-UT  records  of  Duryea  and  Buteo. 

Above,  the  ratios  of  Duryea/Buteo  for  Pg  and  Lg.  Below,  the  ratios 
of  Lg/Pg  for  Duryea  and  Buteo.  The  two  events  have  very  similar 
Lg/Pg  excitation  as  a  function  of  frequency.  Lg  is  dominated  above 
5  Hz  and  only  upper  bounds  are  placed  on  the  Lg/Pg  ratios  above 
5  Hz. 
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Lg/Pg  ratios  are  indicated  above  5  Hz. 

This  comparison  of  these  two  events  indicates  that  although  the  two  events 
differed  in  scaled  depth  by  over  a  factor  of  3,  the  relative  Lg  and  Pg  excitation 
are  virtually  the  same  from  0.25  to  5  Hz.  The  absolute  depths  of  these  events 
differ  by  less  than  a  1/4  the  S-wave  wavelength  at  5  Hz  (about  2000  m/s). 
Because  the  two  hypocenters  were  located  so  close  to  each  other,  the  propaga¬ 
tion  and  station  effects  should  cancel,  and  we  must  conclude  that,  apart  from 
the  spectral  differences  due  to  source  scaling,  the  source  mechanisms  appear 
very  similar. 

In  order  to  better  understand  the  excitation  of  Lg  and  Pg  for  an  explosion 
at  NTS  propagating  in  a  basin  and  range  crustal  structure  we  have  computed 
reflectivity  synthetics  for  some  simple  point  sources.  The  basin  and  range 
structure,  BR3,  is  taken  from  McLaughlin,  et  al.  (1988).  Figure  5.3  shows  the 
response  of  this  structure  to  point  explosions  at  250,  550,  and  660  meters 
depth  as  well  as  a  CLVD  source  at  250  meters  depth.  The  synthetics  are 
broadband  (0  to  2  Hz)  vertical  velocities  from  a  Heavyside  source  function  with 
a  moment  of  10^®  Nt-m.  We  model  spall  by  the  opening  and  closing  of  a  ten¬ 
sion  crack  above  the  explosion  source.  Such  a  tension  crack  source  may  be 
represented  by  the  linear  combination  of  the  moment  tensors  for  a  CLVD  and 
an  explosion.  Consequently,  these  synthetics  can  serve  as  the  fundamental 
elements  needed  to  model  an  axial-symmetric  source  that  includes  both  an 
explosion  and  spall. 
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BR2.  A  =  320  km.  10‘®  Nt-m  CLVD  Source,  h=0.250  km 


0  50  100  150  200  250 

BR2,  A  “  320  km,  10‘®  Nt-m  Explosion  Source,  1  =0.550  km 


Figure  5.3.  Broadband  velocity  (m/s)  synthetics  from  0  to  2  Hz  for  step 
function  sources  with  moment  10^°  Nt-m.  Top  to  bottom,  a  CLVD 
source  at  depth  of  250  m,  an  explosion  source  at  depth  of  250  m, 
an  explosion  source  at  a  depth  of  550  m,  and  an  explosion  source 
at  a  depth  of  660  m. 
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Note  that  the  Pn  and  Pg  amplitudes  are  roughly  the  same  for  the  550  and 
660  meter  explosion  sources  and  that  both  of  these  sources  excite  a  weak  and 
simple  Lg  pulse  that  arrives  at  around  90  seconds.  The  CLVD  source  excites  a 
somewhat  stronger  Pg  and  a  more  complicated  Lg  that  lasts  longer  and  inter¬ 
feres  with  the  early  portion  of  the  Rayleigh  wave.  In  this  layered  seismic  struc¬ 
ture,  the  CLVD  source  excites  Lg  somewhat  more  efficiently  than  the  explosion 
but  also  generates  substantial  Pg.  Therefore,  the  peak  Pg/Lg  ratio  is  not 
greatly  different  for  the  explosion  and  CLVD  source  for  this  structure.  Figure  5.4 
shows  the  theoretical  Lg/Pg  ratios  for  a  vertical  receiver  at  a  distance  of  320  km 
from  three  sources;  an  explosion  at  540  meters  depth,  an  explosion  at  250 
meters  depth,  and  a  CLVD  at  250  meters  depth.  The  Lg/Pg  ratios  were  derived 
by  the  same  procedure  as  those  shown  in  Figure  5.2.  Note  that  one  can  not 
discriminate  between  the  three  sources  based  on  the  Lg/Pg  ratio. 

This  theoretical  result  is  in  stark  contrast  to  the  relative  excitation  of  Lg  by 
CLVD  and  explosion  sources  in  structures  with  high  surface  velocities  (see  Sec¬ 
tion  IV  of  this  report).  Explosions  in  structures  with  high  near  surface  P-wave 
velocities  excite  a  small  amount  of  Lg  when  compared  to  deviatoric  sources 
such  as  a  CLVD  or  tension  crack  source.  CLVD  and  explosion  sources  in 
structures  with  low  near  surface  P-wave  velocities  like  those  at  NTS  produce 
more  nearly  comparable  Pg/Lg  ratios.  It  appears  that  the  additional  shear  wave 
energy  radiated  by  the  CLVD  is  trapped  near  the  surface  by  the  shallow  velocity 
gradients.  These  near  surface  layers  are  more  attenuative  and  the  Lg  is  there¬ 
fore  weaker.  Consequently,  we  conclude  that  the  Duryea-Buteo  test  is  not 
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Lg/Pg  Ratio  of  Green's  Functions 


Frequency  (Hz) 


Figure  5.4.  Lg/Pg  ratios  for  a  vertical  receiver  distance  of  320  km  as  a  function 
of  frequency  for  three  Green’s  functions:  an  explosion  at  550  meter 
depth,  an  explosion  at  250  meter  depth.  The  Lg/Pg  ratios  were 
derived  by  the  same  procedures  used  to  measure  Lg/Pg  and 
Pg/Pg  ratios  in  Figure  5.2.  All  three  sources  have  virtually  the  same 
relative  excitation  of  Lg-to-Pg.  The  frequency  dependence  of  the 
Lg/Pg  ratio  is  largely  due  to  the  differential  attenuation  between  the 
two  phases.  It  is  not  possible  to  discriminate  between  the  three 
sources  based  on  the  Lg/Pg  ratios  alone. 
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sensitive  to  the  differential  excitation  of  Pg/Lg  as  required  to  test  the  hypothesis 
that  spall  is  a  significant  source  of  Lg  at  NTS.  Duryea  may  have  exhibited  more 
spall  than  Buteo  but  theoretically  the  Pg/Lg  ratios  in  the  short-period  band  will 
not  serve  as  a  useful  discriminant  to  test  this  hypothesis. 
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Section  VI 

Research  in  Progress 

In  this  section,  we  describe  two  other  areas  of  research  that  are  in  pro¬ 
gress:  scale  dependent  damage  mechanics  and  finite  difference  modeling  of  the 
Salmon  salt  dome. 

Damage  Mechanics 

Sammis  and  Ashby  (1988)  and  Sammis  (1989)  have  derived  a  model  for 
damage  mechanics  that  is  fundamentally  scale  dependent.  The  strength  of  the 
material  depends  on  the  size  of  cracks  within  the  material  and  the  amount  of 
damage  that  has  been  done  to  the  material.  The  significance  of  this  is  that 
laboratory  measurements  of  rock  strength,  both  measured  directly  and  inferred 
from  small  scale  explosive  tests,  all  show  significantly  higher  rock  strength  than 
that  derived  from  near  field  measurements  of  in  situ  rock.  One  explanation  for 
this  phenomenon  is  that  in  situ  cracks  are  much  larger  than  cracks  in  small 
scale  experiments,  so  in  situ  rock  on  a  scale  appropriate  for  underground 
nuclear  explosions  may  bo  fundamentally  weaker  than  laboratory  scale  rock. 
The  damage  mechanics  model  provides  a  physical  mechanism  with  parameters 
that  scale  to  the  appropriate  size  to  theoretically  simulate  both  laboratory  and 
full  scale  experiments. 

The  basic  relation  for  material  strength  for  this  model  is  Equation  (28)  from 
Sammis  and  Ashby  (1988),  which  gives  the  damage  as  a  function  of  stress 
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State.  This  relation  can  be  reformulated  as  a  function  of  shear  stress,  normal 
stress,  and  damage.  The  function  then  has  the  form: 


(Si -S3)  =  Fi(D)  (Si  +  S3) +  F2(D)  (1) 

where  Si,  S3  are  the  two  normalized  stress  components,  D  is  the  damage,  and 

Fi  and  F2  are  functions  given  by: 
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where  ^  is  the  coefficient  of  friction,  Dq  is  the  initial  damage,  and  a  and  p  are 
empirical  constants  equal  to  approximately  0.7  and  0.45.  The  damage,  D,  is 
defined  by 


D  =  —K(l+aafNv 

w 


(12) 


where  a  is  the  radius  of  the  initial  cracks,  /  is  the  tensile  crack  length,  and 
is  the  crack  density.  The  normalized  stresses  are  defined  by 


S;  = 


cy/(ita) 

K, 


(13) 


where  K"/  is  the  stress  intensity  factor  at  which  the  tensile  cracks  grow  and  a, 
are  the  stress  components. 


Equation  (1)  gives  the  shear  stress  at  which  an  increase  in  crack  growth 
(damage)  occurs.  This  can  be  implemented  within  S-CUBED  finite  difference 
codes  by  checking  at  each  step  in  the  calculation  to  see  if  the  shear  stress  cal¬ 
culated  from  an  elastic  strain  increment  is  below  the  shear  stress  obtained  from 
Equation  (1)  given  the  normal  stress  and  the  previous  damage  in  the  cell.  If  the 
elastic  shear  stress  exceeds  this  value,  then  inelastic  strain  occurs,  causing 
both  an  increase  in  the  amount  of  damage  and  a  decrease  in  the  shear  stress. 
An  additional  constraint  is  required  to  determine  how  much  damage  occurs  on 
each  strain  increment.  This  constraint  can  be  derived  by  relating  the  change  in 
damage  to  the  energy  required  to  increase  the  crack  strain. 
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Simulation  of  Seismic  Radiation  from  the  Tatum  Salt  Dome 

A  second  area  of  continuing  research  is  modeling  the  seismic  radiation 
from  the  Salmon  explosion  using  an  axisymmetric,  linear  calculation  of  an 
explosion  in  a  salt  dome.  The  purpose  of  this  research  is  to  estimate  the 
amount  of  Lg  generated  by  scattering  from  a  complex  earth  structure.  Salmon 
was  well  recorded  and  showed  strong  Lg  signals.  Furthermore,  the  earth  struc¬ 
ture  near  the  source  (Tatum  salt  dome)  is  fairly  well  established.  It  consists  of 
a  high  velocity  salt  diapir  in  low  velocity  sediments  so  it  contains  very  high  velo¬ 
city  contrasts.  This  explosion  therefore  provides  an  ideal  case  study  for  this 
question. 

In  Figure  6.1,  we  show  the  model  that  we  are  using  in  the  calculations  for 
the  salt  dome  structure.  The  compressional  velocity  at  the  source  is  4500 
meters/second.  There  is  a  high  velocity  region  above  the  source  (anhydrite), 
very  low  velocity  sediments  at  the  surface,  and  variable  velocities  gradually 
increasing  with  depth  outside  the  salt  dome.  In  Figure  6.2,  we  show  a  seismo¬ 
gram  from  an  explosion  in  this  structure  at  a  distance  of  16  km.  The  seismo¬ 
gram  was  generated  using  the  two-dimensional  linear,  finite  difference  code 
TRES-2D.  Figure  6.2  compares  the  broadband  (0-3Hz)  seismogram  from  an 
explosion  (depth  850  m)  in  the  2D  structure  with  that  of  an  explosion  in  the  ID 
layered  structure.  The  difference  between  these  2D  and  ID  seismograms  is 
between  10  and  20%  of  the  2D  seismogram  in  this  bandwidth. 


75 


SSS-TR-90-11536 


A  SIMPLE  SALT  DOME  MODEL 

0,0  - - 

I  LOW  VELOCITY  SEDIMENTS 


Figure  6.1 .  Model  used  for  two-dimensional  axisymmetric  modeling  of  the 
Salmon  explosion.  The  P-wave  velocity  at  the  source  is 
4500  m/s.  The  velocity  of  the  anhydrite  layer  is  5800  m/s.  The 
velocity  of  the  sediments  range  from  1 900  m/s  at  the  surface  to 
5000  m/s  at  the  bottom  of  the  grid. 
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2D-1D.  Vertical,  A=16kin  2D-1D.  Radial.  A=16km 


ID  Structure,  Vertical.  A=I61cin  ID  Structure.  Radial.  A=  161cm 


Time  (sec)  Time  (sec) 


Figure  6.2.  Vertical  (left)  and  radial  (right)  0-3Hz  velocity  seismograms  at  a 
distance  of  16  km  for  an  explosion  (depth  850  m)  in  the  2D 
structure  of  Figure  6.1  (bottom)  compared  to  an  explosion  in  a 
layered  structure  (middle).  The  difference  between  the  ID 
layered  structure  and  the  2D  structure  is  shown  at  the  top.  The 
2D- ID  difference  seismogram  is  about  10  to  20%  of  the  2D 
seismogram. 
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